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Abstract. We report on work done in collaboration with Klaus Kinder-Geiger and 
John Ellis which aims at connecting the space-time structure of event generator simu- 
lations with observable output. 



Klaus Geiger was an important driving force in our search for understanding the dynam- 
ics of ultrarelativistic heavy-ion collisions and quark-gluon plasma. An unconventional 
and unique character - this is how we will remember him. Physics was his passion, 
but he also loved his Porsche and his MG (="Ma Geiger"). It is hard for us to accept 
that by his untimely death Klaus, whose research was dedicated so much to the future, 
in particular to the RHIC program at Brookhaven, should be imprisoned forever in the 
past. As his friends and collaborators, we will try to carry on his legacy. 



I INTRODUCTION 

The main focus of Klaus' work during the last years of his life was the event 
generator VNI which describes relativistic heavy-ion collisions and high energy 
particle collisions in terms of a perturbative parton shower Monte Carlo in phase- 
space combined with a simple space-time hadronization prescription. A partonic 
starting point is certainly required at ultrarelativistic collision energies, \fs ^ 100 A 
GeV, where perturbative contributions start to account for a significant part of the 
measured particle production. With his Parton Cascade Model [1], later amended 



by a space-time hadronization algorithm [2] and made publicly available in form of 
the event generator VNI [3] , Klaus Geiger led the way. 

From summer 1997 until Klaus' death we collaborated with him and John Ellis 
on implementing an afterburner for Bose-Einstein correlations into this code and 
testing it in simulations for e + e~ collisions at LEP I and LEP II collisions. Our 
goal was to prepare for an experimental test of the space-time dynamics predicted 
by VNI and to provide theoretical guidance for multi-dimensional Bose-Einstein 
analyses of various event classes with hadronic final states generated at LEP. This 
work [4,5] was interrupted prematurely; our contribution gives an account of the 
present status. 

VNI distinguishes itself from many other high energy event Monte Carlos by 
following the event history in phase-space, not only in momentum space. Space-time 
aspects enter in the numerical simulation in various ways: On the microscopic level, 
the rescattering between produced partons or hadrons is controlled by geometric 
cross sections, and hadronization is modelled geometrically by requiring partons 
to get (in their pair rest frame) closer than 0.8 fm in coordinate space in order 
to form a hadronic cluster. On a macroscopic level this leads to a strong density 
dependence of particle production and absorption rates, further affected by the 
collective expansion of the system which is generated by the rescattering. Finally, 
quantum statistical effects among the produced final state particles, in particular 
Bose-Einstein correlations in momentum space between pairs of identical pions or 
kaons, depend on the phase-space density of the system at the point of decoupling 
and thus on both the momentum-space and space-time structure of the event at 
"freeze-out" (i.e. at the point of the last strong interaction between the particles). 

In this sense all measured quantities in a heavy-ion collision depend to some 
extent on the space-time structure of the reaction zone. The crucial question is, 
however, whether they are sufficiently sensitive to such aspects to allow for a recon- 
struction of the space-time geometry and dynamics. The successful reproduction of 
single-particle yields and spectra by event generators with different space-time fea- 
tures, or none at all like the popular JETSET [6] and PYTHIA [7] generators, seems 
to argue against such a possibility. Two-particle correlations in momentum space, 
however, are sensitive to both the geometric extension of the "fireball" at freeze- 
out and to its collective expansion dynamics [8]. The latter affects the two-particle 
spectra via so-called "x-p-correlations" in the emission function S(x,p). This func- 
tion is the quantum mechanical analogue of the single-particle phase-space density 
of the source at freeze-out: collective expansion correlates the average direction and 
magnitude of the momenta of the emitted particles with their emission points. A 
detailed analysis of 2-pion and 2-kaon correlations in relativistic heavy-ion collisions 
at the Brookhaven AGS and the CERN SPS has recently led to an unambiguous 
demonstration of strong collective dynamics of the fireballs created in these ex- 
periments [9]. On a finer level, however, there remain a number of open physical 
questions whose resolution requires two-particle correlation data of similar quality 
and detail from elementary particle collisions. Such data do not exist, and another 
motivation for our work with Klaus Geiger was therefore to provide stimulation for 



similar experimental analyses of high-statistics high energy data samples like the 
~ 20 million hadronic Z° decays collected at LEP I. 

For sufficiently high secondary particle multiplicities (i.e. at very high energies 
or for large collision systems) it is reasonable to assume that the two particles of a 
selected pair were emitted independently. This allows to express the two-particle 
correlation function in terms of the single-particle Wigner density S(x,p) of the 
source [8]. Neglecting final state interactions (assuming that they can be corrected 
for experimentally [8,9] or theoretically [10] at a later stage), it is given by [8] 

rutr\ . \Sd'xS(x,K)e">'\ 2 \ 
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Here pi,p% are the on-shell momenta of the two particles in the pair while K = 
(pi+p 2 )/2 and q = p\—p 2 are their average and relative 4- momenta. The numerator 
in the second term stands for the product of the measured single-particle spectra 

E p ^ = Jd*xS(x,p) (2) 

where p° = E p = ^rri 1 + pi 1 . The normalization M will be discussed below. 

If, for a given momentum K, the space-time dependence of S(x, K) can be char- 
acterized with reasonable accuracy by a single set of rms widths (i.e. the particle 
emission is not characterized by several widely differing length scales), one can ap- 
proximate S(x, K) by a Gaussian in x. The correlation function then takes the 
simple form [8] 

C(q, K) = N (l + e -9WM*">(*)) , (3) 

where (x^x v )(K) describes the second space-time moments (rms widths or "ho- 
mogeneity lengths") of the effective source of particles with momentum K. Since 
the two measured particles are on-shell, p\ = p\ = m 2 , K ■ q = and only 3 
of the 4 components of q are independent. Different choices for the independent 
components lead to different Gaussian parametrizations of the correlation function 
[8]. We will here use the Cartesian parametrization which eliminates q° = (3 ■ q 
where (3 = K/K° is (approximately) the velocity of the particle pair, q is decom- 
posed into its Cartesian components (q Q , q s , q{) where / denotes the "longitudinal" 
direction (in heavy-ion collisions this is the direction of the beam axis, in e + e~ 
collisions the direction of the thrust axis), o denotes the outward direction, fixed 
by the azimuthal orientation of the transverse pair momentum K± around the l- 
axis, and s denotes the third Cartesian (sideward) direction (defined by K s = 0, 
K=(K ± ,0,K l )). 

Eliminating q° from the exponent in (3) leaves a sum over 6 terms; using fur- 
thermore the azimuthal symmetry of the event sample around the /-axis allows to 
further reduce this to 4 terms involving certain algebraic sums [8] of the rms widths 



(x^x v )(K). So far, however, Bose-Einstein correlations in elementary particle colli- 
sions have been parametrized by much simpler Gaussian forms, involving only one 
or at most two (i^-independent!) size parameters. Such incomplete parametriza- 
tions have the fundamental disadvantage that the fit parameter(s) mix(es) the 
interesting space-time information contained in the rms widths {x^x v )(K) of the 
effective source of particles with momentum K in such a way that they can no 
longer be recovered [11]. An important goal of our work with Klaus was to make 
predictions for the multidimensional shape of the correlation function C(q, K), by 
calculating the complete set of size parameters and their K-dependence. 

The latter is particularly interesting since it signals re- If -correlations in the emis- 
sion function S(x,K). While in heavy- ion collisions the dominant mechanism for 
such correlations appears to be collective expansion of the fireball [11,8,9], contami- 
nations from resonance decays after freeze-out [12] and from x-p-correlations in the 
primary hadron formation process [13] (e.g. from string fragmentation) are known 
to exist. The latter are expected to play a much bigger role in elementary particle 
collisions where a multidimensional analysis of Bose-Einstein correlations may help 
to isolate them. This would provide crucial input into a quantitative discussion of 
Bose-Einstein correlations from heavy-ion collisions where, contrary to elementary 
particle collisions where the multiplicities are much lower, most resonance decays 
can not be reconstructed experimentally and their effects must therefore be simu- 
lated. 

II FROM PHASE-SPACE DENSITIES TO 
MOMENTUM CORRELATIONS 

Event generators like VNI evolve classical probabilities, not quantum mechanical 
amplitudes calculated from properly symmetrized many-particle wave functions. 
The Bose-Einstein correlations among pairs of identical pions or kaons must thus 
be included a posteriori. For notational convenience, we restrict the following 
discussion to a single particle species, say like-sign pions. Let us assume that we 
have generated iVev collision events, and let the mth event (1 < m < N cv ) consist 
of N m such pions in the final state, emitted as free particles from the phase-space 
points {(ri,Pi,ii)}ie[i,N m ]- F° r N ev events, the event generator thus simulates a 
classical phase-space distribution 

1 N cv N m 

Pciass(r, M = F EE ^ (3) (r - fi) 5 (3) (P ~ ft) S(t -U). (4) 

iV ev m= i i= i 

What is needed to calculate two-particle Bose-Einstein correlations according to 
(1) is a prescription which relates this classical phase-space density with the quan- 
tum mechanical single-particle Wigner density S(x,p) of the pion source which is 
supposed to be Monte Carlo simulated by the event generator. In sections II B 
and II A, we focus on two different interpretations of the event generator output 



which we call "classical" and "quantum" , respectively, although these names should 
not distract from the fact that conceptually both are on an equal footing. Each one 
of them leads to a different algorithm for calculating the two-particle correlation 
function from the event generator output: 

{{(^PiMieli,N m] } me[1Nev] =► C(q,K) , (5) 

Before discussing them we list three requirements which any such algorithm should 
fulfill: 

1. Since the event generator evolves classical probabilities, not symmetrized pro- 
duction amplitudes, the generated momenta p iy pj do not show the enhanced 
probability at low relative momenta p { — pj characteristic for Bose-Einstein 
final state symmetrization. The prescriptions (5) should calculate this effect 
a posteriori. 

2. The strength of Bose-Einstein correlations depends on the distance of the 
identical particles in phase-space, not in momentum space. We thus require 
the prescriptions (5) to use the entire phase-space information, and not only 
the generated momentum information. 'Weighting' or 'shifting' prescriptions 
which are based only on the latter [6] may successfully match the measured 
momentum correlations but obviously do not allow to test the simulated space- 
time structure. 

3. In general, Bose-Einstein statistics can affect particle multiplicity distributions 
during the particle production process but classical event generators do not in- 
clude such effects explicitly. Nevertheless they are tuned to reproduce the mea- 
sured average particle multiplicities (N). In order not to destroy this tuning 
we require that the prescription (5) conserves the single-particle multiplicities. 
If the event generator were also tuned to reproduce the average multiplicity 
of identical particle pairs, {N(N — l)/2), i.e. to reproduce not only the mean, 
but also the width of the multiplicity distribution, then the prescription (5) 
should not change that either. Interpreting the correlator as a factor which 
relates the measured two-particle differential cross section to the one simu- 
lated by the event generator, da meas /d 3 pi d 3 p2 = C(q, K) da S i m /d 3 pi d 3 p2, this 
then implies [4,6,14] that Af < 1 in (1). The algorithms discussed below do not 
satisfy this last requirement, i.e. in general they change the width of the multi- 
plicity distribution. However, since the space-time analysis of correlation data 
can be based entirely on the momentum dependence of C(q,K), irrespective 
of its absolute normalization, this does not matter for our purposes. 

A "Classical" interpretation of event generator output 

In the "classical" interpretation [16,5], the phase-space points (ri,Pi,ii) are seen 
as a discrete approximation of the on-shell Wigner phase-space density S(x,p): 



S(x,p) = PclassO,P,£) 



(6) 



This defines both the two-particle correlator via (1) and the one-particle spectrum 
via (2). For a numerical implementation one must replace the momentum-space 
(^-functions in (4) by normalized "bin functions" with finite width e of rectangular 
[16] or Gaussian [5] shape, e.g. 



6^ 



^exp (-(p.- p ) 2 / e 2 ) , 



(7) 



which reduce to the original (^-function in the limit e — > 0. The one-particle spec- 
trum and two-particle correlator then read [16,5] 
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The correlator (9) is the discretized version of the Fourier integral in (1). The 
subtracted terms in the numerator and denominator remove discretization errors 
which would amount to pairs constructed from the same particles. This "classical" 
algorithm, as well as the "quantum" version discussed below, is numerically efficient 
since it involves only 0(N m ) manipulations per event. The calculated observables, 
while being discrete functions of the input, are continuous functions of the measured 
momenta, i.e. no binning in q, K is necessary. 



B "Quantum" interpretation of event generator output 

The "quantum" interpretation [16,15,14,5] of the event generator output starts 
from the observation that, for a given event, i.e. a single term in the sum of (4), 
the simultaneous sharp definition of the particle momenta and positions at emission 
violates the uncertainty relation. In the limit iV ev — * oo the sum is still expected 
to be a smooth phase-space function and, to the extent that the event generator 
provides an accurate simulation of the underlying QCD quantum dynamics, it is 
also expected to respect the uncertainty constraints on any allowed Wigner density. 
Since in practice, however, one has to work with finite numbers of events, one may 
wish to ensure consistency with the uncertainty principle on the single-event level. 

This is achieved [15-17] by associating the centers of single-particle wave packets 
with the set of generated phase-space points {{(ri,Pi,ii)he[i,N m ]}me[i,N ev ]- 
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(10) 



The fi describe quantum mechanically best-localized states, i.e. they saturate the 
Heisenberg uncertainty relation with Axi = a/\/2 and Api = l/y/2a. 

Taking only two-particle symmetrized contributions into account ( "pair approx- 
imation" [14]), all spectra can be written [16,15] in terms of the single-particle 
spectrum Sj(p) corresponding to an individual wave packet at phase-space point i: 
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(12) 




Si (p) = — e" ff ^ . (13) 



Again, the subtracted terms in the numerator and denominator of C(q, K) are finite 
multiplicity corrections which become negligible for large particle multiplicities [15]. 
This algorithm is consistent with an emission function S(x, K) which is obtained 
by folding the classical distribution p c i ass of wave packet centers with the Wigner 
density sq(x, K) of a single wave packet: 

S(x,K) = J ^ficPpidiip^^PiJ^so^x-r^t-i^K-pi), (14) 
s»{x,K) = - z 5{t)e-^ x2 -° 2K2 . (15) 

7T° 

The latter saturates the uncertainty relation with a spatial uncertainty o and a 
momentum uncertainty 1/(7, and the folding ensures that now S(x,K) is always 
quantum mechanically consistent. However, in this algorithm both the one- and 
two-particle spectra depend on the wave packet width a. The role of this parameter 
will be discussed in the context of our toy model study in section III. 

The "classical" and "quantum" algorithms differ only with respect to two points: 

1. The "classical" algorithm has no analogue for the Gaussian prefactor 
exp |o" 2 q 2 /^j in (12) which is a genuine quantum effect stemming from the 
quantum mechanical localization properties of (10). 

2. The Gaussian single-particle distributions Si(p) in the "quantum" algorithm 
are the quantum analogues of the "bin functions" in the "classical" agorithm. 
With the Gaussian bin functions (7) the two agree for the choice a = 1/e. 
Finite event statistics puts a lower practical limit on e in the "classical" al- 
gorithm, but to get accurate spectra one should try to choose e as small as 
possible, by simulating sufficiently many events. In contrast, a in the "quan- 
tum" algorithm is the finite spatial width of the single-particle wave packets, 
and the limit a — > oo (which corresponds to e — > 0) is not physically relevant: 
according to (14) it amounts to an emission function with infinite spatial ex- 
tension and thus to a correlator [15] lim^oo C(q, K) = 1 + S qt Q. 



Ill THE ZAJC MODEL 



Before describing realistic event generator simulations we discuss some analytical 
results for a classical toy model distribution p c i aS s- These are then used to test the 
algorithms of sections II B and II A, by applying those to sets of phase-space points 
{(f i, pi, ^i)}ie[i,Ar m ] generated from the model distribution with a Monte Carlo pro- 
cedure. In this way we can make quantitative statements about i) the a-dependence 
of the "quantum" algorithm, ii) the e-dependence of the "classical" algorithm (es- 
pecially: how small e has to be chosen to extract the limit e — > numerically) and 
iii) the statistical requirements for the algorithms to work. 

The toy model, first introduced by Zajc [18], reads: 

&(r,p,t) = K S(t) exp {-^5y " + J , (16) 

M -= E *w£w *.=wr^. (it) 

The distribution is normalized to a total event multiplicity N. The parameter 
s smoothly interpolates between completely position-momentum correlated and 
un correlated sources. For s — > 0, the position- momentum correlation vanishes and 
we are left with the product of two Gaussians in r and p. In the limit s — > 1 

Jm^frp,*) ~ 5 (3) - -|) S(t) , (18) 

the position-momentum correlation is perfect. The total phase-space volume of the 
distribution V p . s . = (2R S P ) 3 vanishes for s — > 1. This strong s-dependence allows 
to study the performance of our numerical algorithms for different phase space 
volumes. In the following subsections, we discuss the s-dependence of the one- 
particle spectrum and two-particle correlator, focussing in sections III A and III B 
on analytical results, and comparing these in section III C to numerical calculations. 



A The Zajc model in the "classical" algorithm 

Inserting (16) into (6)-(9) one finds [5] 

dN r, N ( P 2 \ 

C(q,K) = l + exp(-qX 2 iassW), (20) 

^ciass(e) = I - e 2/ S ( 2 p2) ( X + 2P 2 (1 - s 2 ) " (2R S P ) 2 ) ' (21) 

We recover the physical HBT radius i? 2 lass = R 2 S {1 - 1/{2R S P ) 2 ) from (21) in the 
limit e — > or by inserting (6) directly into (1). The remarkable fact is that for 
sufficiently large s [18], 



the HBT radius -R^iass becomes negative and the two-particle correlator shows an 
unphysical rise of the correlation function with increasing q 2 . The change of sign in 
(20) seems to be related to the violation of the uncertainty relation by the emission 
function for large s when r and p become strongly correlated. At the critical value 
s crit the phase-space volume V p s , of the source drops below 1. Only for s < s crit 
the distribution p^iass * s a quantum mechanically allowed emission function, i.e. a 
Wigner density. 

The practical importance of the Zajc model in the unphysical limit s > s cr j t is that 
it provides an extreme testing ground for our numerical algorithms. Analytically, 
we conclude already from (21) that in order to be close to the physical limit the 
bin width e has to be small on the scale of the width of the generated momentum 
distribution, 

e^V2P . (23) 

In the above toy model this requirement is independent of the strength of position- 
momentum correlations in the source. 



B The Zajc model in the "quantum" algorithm 

Inserting the model distribution p^^ s into (14) (instead of (6)), we find 
n dN „ N ( p 2 \ 

E ^ P = Ep j2^W 2eXP [-2P 2 )^ (24) 
C(q,K) = l + ex P (-q 2 R 2 qm ), (25) 

R2 ^ R2 ( l -Wn 2 )' w 

In this case, R and P satisfy 2RP > 1 independent of the value of a and, in contrast 
to the "classical" algorithm, the radius parameter R 2 m is now always positive, irre- 
spective of the value of s. Even if the classical distribution p^iJ^ s (f,p,i) is sharply 
localized in phase-space, its folding with minimum uncertainty wave packets leads 
to a quantum mechanically allowed emission function S(x,p) and to a correlator 
with falls off with increasing q 2 as expected. However, the spread of the one-particle 
momentum spectrum (24) receives an additional contribution 1/a 2 . Choosing a too 
small increases this term beyond the phenomenologically reasonable values, choos- 
ing it too large widens the corresponding HBT radius parameters significantly. It 
was argued [15] that a can be interpreted as quantum mechanical "size" of the 



particle, a ~ 1 fm. Given the heuristic nature of these arguments and the signifi- 
cant modifications this implies for the spectra (24) and (25), it is however fair to 
say that presently a mainly plays the role of a regulator of unwanted violations of 
the quantum mechanical uncertainty principle while a deeper understanding of its 
origin in the particle production dynamics is still missing. 



C Numerical simulations in the Zajc model 

We have studied the performance of our Bose-Einstein algorithms by generating 
with a random number generator sets {(^i,Pj, £i)}ie[i,jv m ] of phase-space points 
according to the distribution p^ass an< ^ comparing the numerical results of our 
algorithms to the analytical expressions of section III A and III B. 
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FIGURE 1. Generic properties of the one-dimensional Zajc model, (a): HBT-radius parameter 
(21) of the "classical" interpretation as a function of the position- momentum correlation s. The 
plot shows the HBT radius for different combinations of the model parameters Rq and Po, and 
for different sizes of the bin width e used in the numerical implementation, (b): Same as (a) 
for the "quantum" version (26) of the model. The dependence of the HBT radius on the wave 
packet width a is clearly seen, (c) and (d): The two-particle correlator in the "classical" (c) 
and "quantum" (d) version of the model for different sets of model parameters. The numerical 
results are obtained for 50 events of multiplicity 1000, and agree very well with the analytical 
calculations. 



Fig. 1(a) shows the e-dependence of the HBT radius parameter (21). For fixed 
bin width e, the approximation of the true HBT radius parameter by i?ciass( e ) i s seen 
to become better with increasing P , in agreement with (23). For the HBT radius 
obtained from the "quantum" version of the Zajc model and depicted in Fig. 1(b), 
the situation is both qualitatively and quantitatively different. Now, the HBT 
radius is always positive, since the Gaussian wave packets take quantum mechanical 



localization properties automatically into account. Also, the s-dependence of the 
radius is seen to be much weaker since the wave packets smear out the unphysi- 
cally strong position-momentum correlations present in the classical distribution 
Pciass- The HBT radius depends not only on the geometrical size Rq, and on the 
momentum width P of the source, but also on the wave packet width a. As seen 
in Fig. 1(b), this wave packet width affects the HBT radius and its P -dependence 
significantly for a > Rq. 

In Figs. l(c,d) we present the two-particle correlation functions corresponding to 
these HBT radius parameters for characteristic values of the model parameters. The 
analytical results, obtained by plotting (20) and (25), are compared to the results 
from the event generator algorithms (12) and (9). The plot was obtained using 
50 events of multiplicity 1000. The differences between analytical and numerical 
results originate from statistical fluctuations and become smaller with increasing 
number of events N cv or event multiplicity N m . The good agreement between 
analytical and numerical results in Figs. l(c,d) indicates the relatively low statistical 
requirements of our algorithms. The reason is that both algorithms associate with 
the discrete set of generated momenta p { continuous functions of the measured 
momenta p 1 , p 2 - This smoothens any statistical fluctuation significantly. For the 
"classical" algorithm, small deviations between numerical and analytical results 
are still visible in Fig. 1(c), while the results of the "quantum" algorithm coincide 
within line width. This can be traced back to the Gaussian prefactor exp (— a 2 q 2 /2) 
in (12) which provides an additional smoothening of statistical fluctuations not 
present in the "classical" algorithm. 

IV TWO-PARTICLE CORRELATIONS FROM VNI 

We have applied the "quantum" algorithm discussed above to simulated e + e~ 
collisions at LEP I [19] and LEP II [20] energies from VNI. We have focussed on the 
channels (see Fig. 2) 

e + e~ -> Z° -> qq -> hadrons at y/s = 91.5 GeV , (28) 
e + e~ -> W + W~ -> qq'q'q -> hadrons at y/s = 183 GeV, (29) 

which provide the "cleanest" environment for the study of Bose-Einstein correla- 
tions in high-energy particle collisions. Especially for the first channel there exists 
an impressively extensive and accurate data sample of several million events. For 
high-energy e + e~ collisions, theoretical studies of Bose-Einstein enhancements have 
mainly been performed within the context of the string models [21], which have 
been quite successful in reproducing the distributions of identical particle pairs on 
the level of 1-parameter Gaussian fits of the correlator [22]. Our study [5] is based 
on the event generator VNI and aimed at a full inclusion of the space-time structure 
of the events and a multidimensional correlation analysis. 

The "classical" and "quantum" interpretation of the event generator output dis- 
cussed in section II A and II B provide two different prescriptions for the calcula- 



FIGURE 2. Schematics of the two e + e~ event types (28) and (29): The final-state hadron 
distribution in Z° events (left) is due to exclusively 'endogamous' hadronization of the partonic 
offspring from the qq dijet, whereas in W + W~ events (right) there is, in addition, the possibility 
of 'cxogamous' hadron production involving a mating of partons from the two different W + — ► qq' 
and W~ — > q'q dijets. 

tion of two-particle correlations. This ambiguity in the calculational scheme can be 
traced back to the fact that both algorithms are a posteriori remedies for an incom- 
plete quantum mechanical time evolution which does not account properly for the 
quantum mechanical symmetrization of identical iV-particle states. The ambiguity 
has to be removed by a physical consistency requirement. From our discussion 
of the Zajc model in section III, first crude statements about such consistency 
requirements can be made: 

1. The "quantum" interpretation introduces a wave packet width a which must 
be adjusted to data. The measured HBT radius parameters provide an upper 
bound on a, a < -Rhbt- Hence, the "quantum" interpretation can only be 
consistent with experimental data if the wave packet width is not too large. 

2. In order for the "classical" interpretation of the event generator output to be 
consistent with experimental data, the generated output must not be peaked 
too strongly in phase-space. Otherwise, the HBT radius may be unphysically 
small or even show the wrong sign, cf. section III A. 

The first step in a realistic study of two-particle correlations is necessarily to refine 
these crude statements. To what extent do physical observables depend on the 
choice of algorithm? Which choice of a is consistent with experiment? Are there 
phenomenological reasons to prefer one of the algorithms? These workshop pro- 
ceedings reflect the state of our work when Klaus left us. At the time of the airline 
accident we had only completed the calculation with the "quantum" algorithm, 
and only for a single value of the wave packet width a. A comparison of the two 
algorithms remains to be done. 



A Two-particle correlations at vanishing pair momentum 



Fig. 3 shows the correlator C(q, K) — 1 for different q- values and vanishing pair 
momentum K in the cm. frame of the collision, C(q z , q s , q Q , 0) — 1. The widths of 
the correlator in three Cartesian directions are roughly the same. Quantitatively, 
they are roughly given by the inverse of the wave packet width a which suggests 
that they are dominated by the "quantum mechanical smearing" features of the 
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FIGURE 3. The correlation function of same-sign pions for different values of the relative pair 
momentum q for vanishing pair momentum K, C(q z ,q s , q Q , 0) — 1. 



"quantum" Bose-Einstein algorithm used here. A simulation with the "classical" 
algorithm remains to be done. 

One also sees that near q = the correlation functions show characteristic de- 
viations from a Gaussian shape. However, at non-zero values of the orthogonal 
qr-components, the correlators become nicely Gaussian. The non-Gaussian features 
at small q can be traced to decay contributions from long-lived resonances: Fig. 4 
shows the same correlator with and without such contributions. Neglecting the 
pions from long-lived resonance decays, the correlation function becomes Gaussian 
and wider, reflecting a smaller source size. With long-lived resonances included, the 
effective pion source S(x, K) is larger, resulting in a narrower and non-Gaussian 
correlation function. This clearly demonstrates the sensitivity of the correlation 
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FIGURE 4. The correlator C(q, K) — 1 for K = with (solid lines) and without (dashed lines) 
the contributions of pions stemming from long-lived resonances. 



function on the space-time geometry of the source function S(x, K). 



B Pair momentum dependence of the correlation function 

As discussed before, the dependence of the correlation function on the pair mo- 
mentum K indicates x-i^-correlations in the source and can thus be sensitive to 
the space-time dynamics of hadron production. Fig. 5 shows the correlation func- 
tion C (q, K) — 1 of same-sign pions for various values of the pair momentum 
K = (K L ,K±). Clearly, as K increases, the correlation function becomes wider, 
indicating a smaller effective source, qualitatively (although not quantitatively) 
very similar to the corresponding observations in heavy-ion collisions. This is an 
interesting prediction which to our knowledge has not yet been tested in the LEP 
experiments. According to the more detailed analysis presented in [5] the observed 
Independence can be fully attributed to the /^-dependence of resonance decay 
contributions; however, this may be to some extent an artifact of the employed 
"quantum" Bose-Einstein algorithm whose wave packet width a apparently domi- 
nates the widths of the correlation functions shown above. 



V OUTLOOK 

The rather abrupt ending of the above section on physics results gives a sad 
feeling for the gap left by Klaus. In our collaboration, he was the only one who 
actually knew how to run VNI. Our original motivation came from relativistic 
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FIGURE 5. The correlation function of same-sign pions C(q, K) — 1 for various values of the 
pair momentum K = (Kl,K±), where Kl — K z is the direction along the thrust axis and 
]fi = | if j. | = \Jk \ + is the momentum transverse to it. The correlators are plotted against 
one component of the relative momentum, setting the two other components to zero. 



heavy-ion collisions where the space-time geometry and dynamics of the event plays 
a crucial role in understanding basic measurable quantities. Mainly for this reason 
we developed algorithms which allow to calculate identical two-particle correlation 
functions from an arbitrary event generator output. We tested the accuracy and 
statistical requirements of these algorithms in simple toy model studies and applied 
them in a first realistic study to the hadronic channels in e + e~ annihilation at the 
Z°-peak and near the W + -W~ threshold. However, a comparative study of both 
algorithms is still missing and our main goal, the application of these algorithms 
to the study of event simulations of heavy ion collisions, is not achieved yet. We 
plan to do so in the near future. Also, we plan to make the algorithms described 
in section II publicly available, using the Open Standard for Codes and Routines 
(OSCAR) format, advocated in [23]. 
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